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Abstract 

We study the evolution of the dark energy parameter within the scope 
of a spatially flat and isotropic Friedmann-Robertson- Walker (FRW) model 
filled with barotropic fluid and dark energy. To obtain the deterministic 
solution we choose the scale factor a(t) = vre* which yields a time de- 
pendent deceleration parameter (DP). In doing so we consider the case 
minimally coupled with dark energy to the perfect fluid as well as direct 
interaction with it. 
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1 Introduction 

Observations of distant Supernovae (SNe la) [T] — [5] , fluctuation of cosmic mi- 
crowave background radiation (CMBR) [7J[5], large scale structure (LSS) [51 ITU], 
sloan digital sky survey (SDSS) [HI [12], Wilkinson microwave anisotropy probe 
(WMAP) [13] and Chandra x-ray observatory [TJ] by means of ground and 
altitudinal experiments have shown that our Universe is spatially flat and ex- 
panding with acceleration. This fact can be put in agreement with the theory 
if one assumes that the Universe is basically filled with so-called dark energy. 
The measurement of photometric distances to the cosmological Supernova, sup- 
ported by a number of independent arguments, in particular by the observational 
data on the angular temporal fluctuations of CMBR, shows that the lion share 
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of the energy density of matter belongs to non-baryonic matter. This form of 
matter cannot be detected in laboratory and does not interact with electromag- 
netic radiation. Given the fact that almost three fourth of energy density of the 
Universe originated from dark energy and plays crucial role in the accelerated 
mode of expansion of the Universe, there appear a large number of models ca- 
pable of describing this dark energy. 

Dark energy models with higher derivative terms were constructed by Zhang 
and Liu [TS]. The cosmological evolution of a two-field dilaton model of dark 
energy was investigated by Liang et al. [16] . Viscous dark energy models with 
variable G and A were studied by Arbab [T7] . The modified Chaplygin gas with 
interaction between holographic dark energy and dark matter was investigated 
in Ref. [18] . The tachyon cosmology in interacting and non- interacting cases in 
non-flat FRW Universe was studied in Ref. [15] . In this Letter we study the 
evolution of the dark energy parameter within the framework of a FRW cosmo- 
logical model filled with two fluids. In doing so we consider both interacting 
and non-interacting cases. 



2 The Metric and Field Equations 

We consider the spatial homogeneous and isotropic Friedmann- Robertson- Walker 
(FRW) metric as 



ds 2 = -dt 2 + a 2 (t) 



dr 2 



1 - kr 



(1) 



where a(t) is the scale factor and the curvature constants k are — 1,0, +1 for 
open, flat and close models of the universe respectively. 

The Einstein's field equations (with 8ttG = 1 and c = 1) read as 

IK \M t;. (2) 

where the symbols have their usual meaning and T( is the two fluid energy- 
momentum tensor consisting of dark field and barotropic fluid. 

In a co-moving coordinate system, Einstein's field equations ^ for the line 
element (J]) lead to 



Ptot = - 2- + — + — , (3) 
\ a a gt / 

and 

J a 2 k\ 
, tot ^- + -J, (4) 

where ptot = Pm + Pd and ptot = Pm + Pd- Here p m and p m are pressure and 
energy density of barotropic fluid and po & Pd are pressure and energy density 
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of dark fluid respectively. 

The Bianchi identity G'?- = leads to T'A = which yields 



ptot + 3- (ptot + Ptot) = 0. (5) 



The EoS of the barotropic fluid and dark field are given by 

_ Pm 

W m — , 

Pm 

and 



respectively. 



(6) 



ud = — , (7) 

PD 



In the following sections we deal with two cases, (i) non-interacting two-fluid 
model and (ii) interacting two-fluid model. 



3 Non-interacting two-fluid model 

First, we consider that two-fluids do not interact with each other. Therefore, the 
general form of conservation equation ([5]) leads us to writing the conservation 
equation for the dark and barotropic fluid separately as, 



and 



Pm + 3^ (p m + Pm) = 0, (8) 



pD + 3-(p D +p D ) =0. (9) 
a 



Integration of Eq. ([5]) leads to 

p m =p n a-^ 1+ ^\ (10) 

where po is an integrating constant. By using Eq. (fTUf in Eqs. ^ and ((U), we 
first obtain the po and po in term of scale factor a(t) 



a 2 



and 

= (z- + - 2 + K )- PoU> m a- 3{1+Um) . (12) 

Now we take following ansatz for the scale factor, where increase in term of 
time evolution 

a(i) = Vie*. (13) 
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Figure 1: The plot of EoS parameter vs t for p = 10 and oj m = 0.5 

The motivation to choose such scale factor is behind the fact that the universe 
is accelerated expansion at present and decelerated expansion in the past. Also, 
the transition redshift from deceleration expansion to accelerated expansion is 
about 0.5. Thus, in general, the DP is not a constant but time variable. By the 
above choice of scale factor yields a time dependent DP. 




By using this scale factor in Eqs. (|TT|) and (TT2)) . the pn and pn are obtained as 

~(l + t) 2 k 



Pd = 3 



U 2 



te 1 



poitey 



and 



PD 



3(1 + t) 2 - 4 k 



At 2 



+ 77+ Pow m (te* 



(14) 



(15) 



respectively. By using Eqs. (IT41 and (fTS")) in Eq. (O, we find the equation of 
state of dark field in term of time as 



3(l+t) 2 -4 . fc , /. t \-|(l+ w m ) 



(l+t) 2 , Jl 
4t 2 "f" te* 



- Po(*e*) 2 



n-|(i+^ m ) 



(16) 



The behavior of EoS for DE in term of cosmic time t is shown in Fig.l. It is 
observed that although for open, close and flat universes the EoS parameter 
is an increasing function of time, the rapidity of its growth at the early stage 
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depends on the type the universe. Later on it tends to the same constant value 
independent of the types of the universe. 



The expressions for the matter-energy density fl m and dark-energy density £Id 
are given by 

r> P m a t\—%(i+w m ) tT7\ 



and 



^ Pd , 4fci At 2 , f . 3,-,,., \ , , 

respectively. Equations ([T7)l and (fT8|) reduce to 

4fci 

n = n m + n 1 , = i + 5rrTW . (19) 

From the right hand side of Eq. (TIT)]) it is clear that in fiat universe (fc = 0), 
= 1 and in open universe (k = — 1), O < 1 and in close universe (fc = +1), 
f2 > 1. But at late time we see for all flat, open and close universes il — > 1. This 
result is compatible with the observational results. Since our model predicts a 
flat universe for large times and the present-day universe is very close to fiat, the 
derived model is also compatible with the observational results. The variation 
of density parameter with cosmic time has been shown in Fig. 2. 

We define the deceleration parameter q as usual, i.e. 



aa a .„„. 

9 = -7a=-^m- ( 2 °) 



Using Eqs. §5§ and (gj), we may rewrite Eq. (|2"0j) as 

1 



[p m (l+3u m ) + p D (l + 3u D )}. (21) 



1 6H 2 

On the other hand, using Eq. (fT3j) into Eq. (|20)) . we find 

2 



(1 + *) 2 



1. (22) 



From Eq. ([22]). we observe that q > for t < 0.41 and q < for t > 0.41. 
This behavior of q is clearly depicted in Fig. 3 

A convenient method to describe models close to A CDM is based on the 
cosmic jerk parameter j, a dimensionless third derivative of the scale factor 
with respect to the cosmic time [20] — [23 . A deceleration-to-acceleration tran- 
sition occurs for models with a positive value of jo and negative qo . Flat A 
CDM models have a constant jerk j = 1. The jerk parameter in cosmology is 
defined as the dimensionless third derivative of the scale factor with respect to 
cosmic time 

m - jp\- (23) 
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G 



and in terms of the scale factor to cosmic time 



m - ^ 

where the 'dots' and 'primes' denote derivatives with respect to cosmic time and 
scale factor, respectively. The jerk parameter appears in the fourth term of a 
Taylor expansion of the scale factor around ao , 

^ = 1 + H (t - 1 ) - \q Q Hl{t - 1 ) 2 + \joH*{t - t Q f + 0[(t- t ) 4 ] , (25) 
ao zd 

where the subscript shows the present value. One can rewrite Eq. (|23|) as 

j(t) = q + 2q 2 -jj. (26) 
Equations (|2"21 and (|2"6"|) reduce to 

This value overlaps with the value j ~ 2.16 obtained from the combination of 
three kinematical data sets: the gold sample of type la supernovae [53], the 
SNIa data from the SNLS project [25], and the x-ray galaxy cluster distance 
measurements 126 at t ~ 0.05. 



4 Interacting two fluids model 

Secondly, we consider the interaction between dark and barotropic fluids. For 
this purpose we can write the continuity equations for dark fluid and barotropic 
fluids as 

prn + 3- (p m + Pm) = Q, (28) 

and 

Pd + 3- (pa + Pd) = —Q- (29) 
a 

The quantity Q expresses the interaction between the dark components. Since 
we are interested in an energy transfer from the dark energy to dark matter, we 
consider Q > 0. Q > ensures that the second law of thermodynamics stands 
fulfilled [37] • Here we emphasize that the continuity Eqs. (|28p and ([2U]) imply 
that the interaction term (Q) should be proportional to a quantity with units 
of inverse of time, i.e Q oc j . Therefore, a first and natural candidate can be 
the Hubble factor H multiplied with the energy density. Following Amendola 
et al. [28] and Gou et al. [29], we consider 

Q = 3Ha Pm , (30) 
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where a is a coupling constant. Using Eq. (|50|) in Eq. ([25]) and after integrating 
the resulting equation, we obtain 



Pn 



p a 



(31) 



By using Eq. ([3~Tj) in Eqs. (J3|) and (@J, we again obtain the pn and in terms 
of scale factor a(t). 



Pd = 3 



and 



p a 



-3(H 



-<r) 



-3(l+w m .-cr) 



(32) 



(33) 



respectively. Putting the value of a(t) from Eq. (jTU)) in Eqs. (|52"j) and (|33l) . we 
obtain 

^(1 + t) 2 fc 



PD 



and 



4t 2 te* 
3(1 + <) 2 - 4 fc 



Po{te l 



4t 2 



(34) 



(35) 



respectively. Using Eqs. ([34]) and (1351 in Eq. (O, we can find the EoS parameter 
of dark field as 



3(l+t) 2 -4 , fe , /. tw|(l+u m -v) 
172 r JTt + PoW m \te ) 2\ mi 



(1+*) 2 
4t a 



p (te*)-i( 1+w '"-' J ) 



(36) 



The behavior of EoS in term of cosmic time t is shown in Fig. 4. It is observed 
that like the minimal coupling case, the EoS parameter is an increasing function 
of time for all close, open and flat universes, the rapidity of its increase at the 
early stage depends on the type of universe. At the later stage of evolution it 
tends to the same constant value independent of the types of the Universe. The 
EoS parameter of DE begins in phantom region and tends to —1 (cosmological 
constant. 



The expressions for the matter-energy density VL m and dark-energy density 
d are given by 



and 



n, 



Pd 



Pv 



At 2 



3H 2 3(1 + t) 2 
4kt 



1 



p (te 



At 2 



3H 2 3(l + i) 2 e* 3(1 + t) 2 

respectively. From Eqs. (|37| and (|38|) . we obtain 



poite 1 ) 



t\-3n(l+ui m -cr) 



n, 



i 



3(1 + t)V 



(37) 



(38) 



(39) 
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Figure 4: The plot of EoS parameter vs t for p — 10, ui m — 0.5 and a — 0.3 



which is the same as Eq. (JT9J) . Therefore, we observe that in the interacting 
case the density parameter has the same properties as in the non-interacting 
case. The expressions for deceleration parameter and jerk parameter are also 
the same as in the case of non-interacting case. 



Studying the interaction between the dark energy and ordinary matter will 
open a possibility of detecting the dark energy. It should be pointed out that 
evidence was recently provided by the Abell Cluster A586 in support of the in- 
teraction between dark energy and dark matter [30l [31] . We observe that in the 
non-interacting case only open and flat universes can cross the phantom region 
whereas in interacting case all open, flat and close universes can cross phantom 
region. 



5 Concluding Remarks 

In summary, we have studied the system of two-fluid within the scope of a spa- 
tially flat and isotropic FRW model. The role of two-fluid minimally or directly 
coupled in the evolution of the dark energy parameter has been investigated. In 
doing so the scale factor is taken to be an exponential law function of time. It 
is concluded that in the non-interacting case only open and flat universes cross 
the phantom region whereas in the interacting case all three universes can cross 
the phantom region. 
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